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^ I Abstract: The identification of the correct model for physics beyond the Standard 

^yJ [ Model requires the determination of the spin of new particles. We investigate to 

O I which extent the spin of a new particle X can be identified in scenarios where it decays 

dominantly in three-body decays X — )■ ff'Y. Here we assume that y is a candidate 
for dark matter and escapes direct detection at a high energy collider such as the 
'k>( I LHC. We show that in the case that all intermediate particles are heavy, one can get 

H ' information on the spins of X and Y at the LHC by exploiting the invariant mass 

distribution of the two standard model fermions. We develop a model-independent 
strategy to determine the spins without prior knowledge of the unknown couplings 
and test it in a series of Monte Carlo studies. 
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1. Introduction 

With the start of the Large Hadron Collider (LHC) the direct exploration of physics 
at the TeV scale has begun. The hunt for new physics is one of the major topics 
in the experimental program of the LHC. Many of these models predict partners of 
the known standard model (SM) particles, which usually have the same quantum 
numbers and properties but for the mass and the spin assignment. For example, in 



supersyminetric (SUSY) models the fermions have scalar partners whereas in models 
with universal extra dimensions (UED) fermionic partners are predicted. Due to the 
astrophysical requirement of explaining the dark matter (DM) relic density of the 
universe, these models usually invoke an additional discrete symmetry leading to a 
new stable particle which in general escapes detection at future collider experiments. 
Examples are R-parity in SUSY models or Kaluza-Klein-parity in UED where e.g. the 
lightest neutralino or the lightest Kaluza-Klein (KK) excitation of the vector bosons 
is the corresponding DM candidate, respectively. The generic signature at LHC 
are in both cases SM-fermions with high transverse momentum and missing energy 
stemming from the escaping DM candidate. 

An important question is: How can one distinguish between different models? These 
models differ in the spins of the predicted new particles and, thus, one has to develop 
methods to get information on the spin. First attempts have been made for s-channel 
resonances |1], 0, ^ ^ and in case of subsequent two-body decays of the new particles 
[|, §, 0, g |, |TD|, ITT], 0, ig, 0, [T5|, ig [T^] where in many cases model dependent 



assumptions had been made. An additional possibility to get information on the 
spin is cross section measurements provided one knows the representation of the 



particle produced |jT8[, e.g. whether it is a colour triplet or a colour octet. Hardly 
any attempt has been made so far in case of three-body decays but for the case 
of distinguishing a gluino from the KK excitation of a gluon |T^ and the quantum 



interference method |^0|, |21|, g^, ^ which, in principle, is also valid for three-body 
decays. The quantum interference method requires fully reconstructed events which 
can be achieved at the ILC and only in few processes at the LHC. In contrast, our 
method in this paper does not require full reconstruction of events. In this paper we 
start a series of investigations on how one can extract information on the spins of 
new particles in a model independent way if three-body decays are dominating. 
We concentrate here on the case of the direct production of a new particle X de- 
caying via a three-body decay into two SM-fermions and a new invisible particle F, 
which escapes detection. We will show that the invariant mass distribution of the 
two detectable fermions contains sufficient information to extract the spins of the 
unknown particles in such a decay. 

In this paper, we consider cases where the intermediate particles are very heavy 
compared to the decaying one. Examples of such cases are e.g. split SUSY with very 



heavy scalars [^, split UED [^ or higgsless supersymmetric models [^. As it 



turns out, in this limit it is possible to determine the spin of the decaying particle 
and the invisible particle provided the masses of the new particles are known to some 
extent. In contrast to model dependent approaches to spin determination |Jl9|, ^, we 
follow a bottom-up approach that can be used for all models with heavy intermediate 
particles and only renormalizable operators. The case of lighter intermediate particles 
will be presented in a subsequent paper. 
This paper is organised as follows: We first present the basic ideas using a toy 



model in section |^. Afterwards, we present a convenient parametrisation of the 
differential width for the three-body decays as the product of a phase space factor and 
a polynomial. In section | we develop a strategy to determine the spins exploiting 
specific properties of the coefficients of this polynomial. We also investigate the 
impact of different decay topologies and discuss qualitatively the influence of the 
mass of the intermediate particle. In section ^ we test our strategy with the help 
of Monte Carlo examples and conclude in section |^. The formulas of the various 
coefficients are given in the appendix. 



2. Basic Idea and General Setting 

We investigate decays of the type X — )■ f fY where X and Y are new particles being 
either scalars, vector bosons or fermions. Here we assume that all 2-body decays 
of X are either kinematically forbidden or at least loop-suppressed compared to the 
tree-level three-body decays considered. As mentioned above we assume that all off- 
shell particles, which we denote collectively by /, mediating these decays are much 
heavier than X, e.g. mi ^ mx- In practice it is sufficient to assume mi > 5mx as 
we will show below. In addition we assume that F is a colour singlet as it should 
serve as potential dark matter candidate. 

We will be as general as possible by taking the most generic Lagrangian with arbitrary 
couplings of 0{1) and dimension 4 operators. From this we calculate the widths for 
the decays of Y assigning different spins to X and Y, respectively. To simplify the 
notation we will abbreviate the decays S — )■ f fS, S — )■ ffV, V — )■ ffS, V — )■ ffV 
and F -^ ffF by (5, S), {S, V), {V, S), {V, V) and (F, F), where S, V and F stand 
for scalar, vector boson and fermion, respectively. Note, that the fermionic case 
covers both, Dirac- and Majorana-fermions. 

After integrating over the momentum of the escaping particle Y, we expand the 
differential widths in powers of e = m,x/iTii and give the resulting expressions as a 
phase space factor times a Laurent series (actually polynomials in most cases) of a 
dimensionless quantity s which is derived from the invariant mass s = {pf -|-pj)^, pf 
and Pf are the four momenta of the SM- fermions. Note, that s = m?r7 in the case of 
subsequent two body decays studied in the literature [^, ^ 0, |], |], |1^, ^, |2^ . 
Most of the features can be understood by considering the decays of a particle X 
charged under SU{?>) into two massless quarks and Y . Considering coloured particles 
in the first place is motivated by the fact that they in general have sizable cross 
sections at the LHC. Moreover, due to gauge invariance only a subset of all topologies 
are allowed which simplifies the obtained expressions considerably. Therefore we will 
first discuss these cases. The additional features of either taking X as an SU{?>) 
singlet and/or the case that the SM-fermions being massive (i.e. top quarks) will be 
discussed afterwards. 



The coefficients in the Laurent series depend obviously on the couphngs and masses 
of the particles involved and one might ask if and how one gets information on 
the spins without knowing these quantities. It turns out we have to assume in our 
approach that the mass of Y and the mass difference {mx —iTiy) are known within a 
given uncertainty but in general we do not need any information on the underlying 
couplings. The basic idea is that different spin assignments lead to different relations 
between these coefficients which can be exploited. There is however one obstacle: one 
cannot exclude on logical grounds that there is a 'conspiracy' between the couplings 
suppressing the dominant terms in the e expansion. This complicates life somewhat 
but even in that case relations between the coefficients are maintained as discussed 
below. 

2.1 Dependence of the invariant mass distribution on the spins of un- 
known particles 

We ffist discuss a set of toy models where the new particles X and Y are either scalars 
and/or vector bosons coupling to massless w-quarks and an additional heavy fermion 
which we assume to be a Majorana-fermion. The invariant fermion mass squared is 
s = {pf + 'PfY = E1E2 ■ {1 — co?,6fj) where Ei are the SM fermion energies and Ofj is 
the angle between them in the rest-frame of X. In Figure we show the differential 
decay width divided by the phase space factor which is shown independently (red/full 
line). The behaviour of the different curves can be understood using helicity and spin 
arguments. 

The SM-fermions have definite helicity states as we have assumed them to be mass- 
less. There are two kinematical configurations where the spatial angular momentum 
in the rest-frame of X is zero corresponding to cosOj^j = ±1: (i) The particle Y 
is at rest and the two fermions are back to back corresponding to cos6'jj = — 1 
with s = Smax- Ih this case the total spin of the fermions sums up to one if it is 
either an ulUl or urUr combination whereas the total spin of the fermions is zero 
for the ulUr and urUl combinations. Here we have introduced ul,r = Pl,r u with 
Pl,r = (1 -F 75)72. (ii) The particle Y flies opposite to the two fermions which are 
flying parallel now corresponding to cos 9ff = 1 with s^m = 0. In this case the total 
spin of the fermions sums up to one if it is either a ulUr or a urUl combination 
whereas the total spin of the fermions is zero for the ulUl and urUr combinations. 

iS,S): 
The matrix element Mfi for this decay has the generic form 

Mfi = Uu {grPR + giPL)Si{n*iPR + nlPL)vu ^ — Uu{grn*iPR + gin*.PL)vu 

corresponding to the ulUr and urUl combinations since in the limit of a very heavy 
intermediate particle the corresponding propagator reduces to Sj = l/{^ — rrij) ~ 



as 




Figure 1: Differential width divided by a pliase space factor PS for the different decays 
X -^ ffY, X,Y G S,V taking mf = 0, TUx/my = 0.1 and all couplings equal. In addition 
the phase space factor is drawn. 

l/mj. As the total angular momentum of the final state has to be zero, the matrix 
element has to vanish in case cos Off = 1 (s = 0) as can also be seen in Fig. ^ whereas 
for cosOff = —1 the helicity assignment yields a non- vanishing matrix element. In 
the plot we have taken for all cases rii = 71^. = gi = Qr- 

{s,vy. 

This process has a more involved structure since the vector boson can have polari- 
sation ±1, 0. Using the matrix element one sees that in contrast to the previous case 
one expects a non- vanishing matrix element for all s 

Mfi = Uu {grPa + giPL)Sn^{n*iPR + n*PL)vueY 
- — Uu^iiigrKPi + gin*iPR)vueY- 

Till 

These are the ulUl and urUr combinations where for s = the total spin of the 
fermions is zero and for s = Smax it is one. In the limit Of-j — ;■ (s — > 0) the amplitude 
is proportional to {{mx / frtyY ~ 1) {fn\ — fTtY)- The first factor is larger and diverges 
as ruY — ;• which reflects the longitudinal component of the vector boson. This can 
be nicely seen in the blue (dashed) line in Fig. |I|. 

[v^sy. 

In this case, we start with a spin 1 boson. A similar reasoning as before shows 
that the general trend should be opposite to the S — )■ f fV case which is confirmed 
by Fig. 1^ (pink/small dashed line). 

iy,vy. 

This decay can have several helicity combinations which we want to sketch here. 
The matrix reads 



Mfi = Uu l^l{,grPR + giPL)Snu{n*PL + nlPR)VueYe 



X 



- — Uu l^luigrKPL + gin'^ Pi)vueYex 

As in the (5*, 5*) case we find the ulUr and urUl combinations. However, now in 
principle for all s one can expect a non-vanishing matrix element squared. However, 
as can be seen in Fig. |1| it can be zero for s = for special combinations of the 
couplings, e.g. as in our case with rii = n^ = gi = gr (cyan/dot-dashed line). In this 
sense the example shown is an extreme case and in general the matrix element will 
be non-zero for s = 0. 

2.2 General structure of the differential widths 

We now discuss the general structure for the decays 

X{p,mx) -^ f{qi,mf) +7(^2,"^/) + l^(g3,"^y) (2.1) 

with X and Y being scalars (S), vectors {V) or fermions (F) and heavy interme- 
diate particles / with mass mj. We will consider several contributions at the same 
time and assume that the masses of all intermediate particles are equal to maximise 
interference effects which usually complicate things. Beside the usual Mandelstam 
variables 

s = {p- qsY = (gi + g2)^ ; t = {p- ga)^ = (gi + gsf 

and u = {p — qiY = (g2 + QsY = —s — t — m\ — my + '2.m?f 
we introduce the dimensionless parameters Tj and s 

. (4r| + {ry - 1)^) - ^ ^^ ^^ Mc ,, ,, 

where Mc denotes dimensionful couplings, e.g. as they appear in the ZZH vertex 
or the trilinear soft SUSY breaking parameters Ai. We note that Smin = ~1 and 
s =1 

■J max -'- ■ 

We expand the matrix elements squared in powers of e = mx/Tni and perform an 
integration over t as in this way we integrate over the momentum of the unobserved 
particle Y . For the t-integration we find the boundaries 

t± = \rn\ ( (ry + lf + s {4r] - {ry - if) 



± ((r I? Ar-\ X a-^^)(a-^)((^^-l)^-4r|) + 8ry) \ 



PS = f^dt 



where we have also defined the 'phase space' PS. In this way the differential decay 
rate can be written as 
dV PS ( Z A „ „ „ „ „. z. .3 , .. .4 



+ — — - ^B^C -s^D-r + E-s^ + F-s 



ds (27r)3 256 mx \ {as + hf as + b 

(2.5) 

where a = ({ry - if - 4r|) and b = [{ry - 1)^ + 4r|) . The prefactors Z,A,...,F 
are functions of e, the Tj and the couplings. Note that Z and F only appear in the case 
of the decay {V,V). As we exemplify the main features for massless SM-fermions, 
we display Eq. (pl^) for this case: 



t± = Wx {i-ry + If - s{Ty - if ± (1 - Ty) ^(1 " S) ((1 - S) (ry - 1)2 + 8ry) 

PS = ImUl - Ty)^il-s){{l-s){Ty-lf + 8Ty) 

(2.6) 

2.2.1 Decays of bosons 

As one can see from Table 0, there are three 'topologies' which contribute differently 
to the decay rate. The second 'topology' only contributes if both, the Y- and the 
X-particle, are their own anti-particles. Obviously, topologies 1 and 2 will in general 
contribute at ^(e^) whereas the third one in general only at 0{e^) due to the different 
structures of the propagators. Only in the case where the dimensionfull scalar- vector- 
vector or triple scalar couplings are of 0(171 j), the third topology might contribute at 
a smaller power of e as will be discussed below. In the further calculation we neglect 
terms higher than 0{e^). 
The generic Lagrangian density for these decays reads as 

A,,,, = XjfGif + YjfNJ + Ijfif + hX.Ykfijk + h.c. (2.7) 

where i,j, k = s,v and Jj is the intermediate off-shell particle. The generic couplings 
are given by 

dr.Gs = {g{r, s)Pr + g{l, s)Pl) ; G, = r ioir, v)Pr + g{l, v)Pl) 

Nr. N, = {n{r, s)Pn + nil, s)Pl) ; N, = ^ {n{r, v)Pr + nil, v)Pl) (2.8) 

Tr. Ts = {s{r)PR + s{1)Pl) ; T, = ^ {v{r)PR + v{1)Pl) 

and for Vijk 

r,,, = c{s)Mc] r,,, = r,,„ = ^(^;^ v'^, s) = c{v)Mcg^'' 

Vs., = Tis,v^,v'') = c{s)Mcg^'^- T.^.s, = T{v^',si,S2) = c{v){ps, - Vs,Y 

Ts^S2V = ^s,vs2 = T{si,S2,V^') = C{s){ps2 - PsrY 

r„,,,,3 = V{vl, vl v^^) = c{v){{p,, - p,,rg''P + (p., - p^^fg^P + {p,, - p.^Yg^ 

(2.9) 
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Table 1: Topologies for the decays of bosons X — t- f fY (top) and fermions (F,F) (bottom) 
with i, j G {l,r} (see also Eq. (p^ ). 



where the indices of the vertex expressions Vijk correspond to those of the Lagrangian 
in Eq. O. 



The matrix element including all topologies shown in Table |T] reads 
Mij = u{qi,mf)Gi^p 7°A^J7° v{q2,mf)ei{p)ej{q3) 
+u{qi,mf)Nj:^p 7°Gj7° v{q2,mf)ei{p)ej{q3) 
+ Yl ^kijWp^iu{qi,mf)-f^Tl-^^v{q2,mf) 

k=s,v 

with the 'polarisation' vectors 



e^ = 1; e^=?' 



(2.10) 



and the fermion and boson propagators 






P~ 'p^-mf ' 



^'* p'^—m'j ' 



~5 2~ 



(2.11) 



2.2.2 Decays of fermions 

In this case the generic Lagrangian is given by 

£, = hjG^M,, + lijNiMy + hWyVM. + hlT.f + h.c. 



(2.12) 



where M^^y denote the spinors of the new fermions and i = s,v denotes whether the 
exchanged particle is a scalar or a vector boson. The couplings are similar to those 
above: 



M,, 



+ 
+ 



u{qi,mf)Giu{p,mx) 
u{qi,mf)Niv{q3,mY 
u{qs,mY)Tiu{p,mx) 



W, 



p,i 



w. 



p,l 



w. 



p,l 



u{qs, my)7°A/'/7S(g2, "^/) 
v{p,mx)'J°Gl'j°v{q2,mf) (2.13) 

u{qi, m/)7°7;^7°f (g2, ruf) 



corresponding to the topologies given in Table ^ with the same couplings as in Eq. 



!8|, 2^) and additionally: 



i = s: T, = idis,r)PR + d{s,l)PL) 
i = v: T, = Y{d{v,r)Pn + div,l)PL) 



(2.14) 



3. Strategy for Spin Identification 



In this section we discuss the strategy for discriminating the various scenarios with 
different spins assigned to the particles X and Y. The procedure is to find suitable 
relations or to mark the signs of the different coefficients Z,A,...,F. This is done 
in the second part of this chapter. Before this we will have a look at the different 
topologies and their contribution depending on the chosen colour structure. 



Our main focus here is on the case of massless fermions as this is aheady sufficient to 
get the required information. This immediately imphes a considerable simplification 
because some of the coefficients for the differential width are zero and we obtain 

^ = 7T^. {B + Cs + Ds^ + Es' + Fs') (3.1) 

as [Zn)'^ 256 nix 

The formulas for the coefficients are given in the appendix for the corresponding 
lowest order in e. For the case {S, S) we also give the higher orders up to e^. More- 
over, it turns out that the decays in top-quarks, the only SM fermion with mass of 
(9(100 GeV), behave in the same way and the discrimination is also possible in this 
case as will be discussed at the end of this section. This implies that one has a second 
system to test the spin assignments in an independent way. 

We start with a subset of the topologies given in Table |l], namely topologies 1 and 
2 as these are typically realised in extensions of the SM, e.g. in models with extra 
dimensions or in SUSY. Moreover, we will first further restrict ourselves to scenarios 
where X is charged under SU{3), e.g. a colour octet gluino or a KK excitation of a 
gluon, and Y is electrically neutral and uncharged under SU{3). This is motivated 
by the fact that the LHC is a hadron collider. In the second step we add the third 
topology. But it turns out that in the case of scalar contributions to {S, S) and 
{S, V) gauge invariance allows only two additional terms because the SM-fermions 
are charged under SU{3). However, these contributions will in general be of order e'^ 
due to the boson propagator except for the case where the trilinear scalar coupling is 
of order mi in which there might be contributions at order O(e^). In the third step 
we will also discuss the complications and their potential solutions in case that X is 
an SU{3) singlet. 

3.1 Signs of the coefRcients 

It turns out that some of the coefficients have a definite sign independent of the 
couplings and masses involved. This important fact will be used later to discriminate 
between the different spin assignments of X and Y. We collected the signs of the 
different coefficients for all decays of bosons in Table |^ where we have expanded the 
coefficients in powers of e, e.g. 

4 

B = J2Bke' (3.2) 

k=2 

Some of the signs in Table ^ are obtained analytically but several are gained numer- 
ically by scanning and inserting random couplings in the range [—1,1]. There are 
some coefficients where the sign cannot be determined without knowing the mass 
ratios or the couplings which are marked by " ±" . Moreover, we have put a when- 
ever the coefficient itself vanishes. We give the signs for three cases, ordered from 
the most general one to the most restricted one: (i) X is an SU{3) singlet, where 
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Table 2: Signs of the coefficients for the case of a boson decaying into another boson and 
massless SM-fermions in the final state for different powers of e. The rows correspond to 
the cases: (s) X is an SU{3) singlet, (c) X is charged under SU{3) taking all possible 
topologies into account and (1 + 2) X is charged under SU{3) taking topologies 1+2 into 
account. The ± marks the cases where the sign cannot be determined without knowing 
the masses/couplings and marks the cases with a vanishing coefficient. 

all topologies of Table |l] contribute. The corresponding columns are denoted by s. 
(ii) X is charged under SU{3) and all possible contributions allowed are taken into 
account and the corresponding columns are denoted by c. (iii) X is charged under 
SU{3) and only topologies 1+2 contribute and, thus, the corresponding columns are 
denoted by 1+2. In the subsequent sections these cases will be discussed in the re- 
versed order focussing on the terms of order e^. The e^ order is only of interest, if 
the leading order is zero, which is the case for some special coupling arrangements 
discussed in section R^. 



In case of a new fermion X it turns out that the result does neither depend on the 
spin of the exchanged particle nor on the topology, e.g. it does not matter if all 
topologies are taken or only a subset. Since we have only bosonic propagators there 
are only the 0(e'^) contributions and we find: 



sign(54 



+ 



sign(C4 



± 



sign(D4 



sign(E4 







(3.3) 



3.2 Decays into massless SM-fermions in case of topologies 1+2 

Let's assume that we have measured the differential decay width of a new particle 
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and determined the coefficients introduced above accurately in a ffi. In section ^ 
we will discuss first Monte Carlo studies at the parton level where we also review 
the obtainable accuracy. This can be combined with our knowledge on the various 
coefficients introduced so far to determine the spins of the new particles or at least 
to exclude certain possibilities. The main strategy is summarised in Fig. ^ and 
explained in some more detail below. 

Let's start with the E term which is only non-zero in the (V, V) case. This imme- 
diately implies that (V, V) is preferred once the 'measured' E term is larger than 0. 
For consistency we check that B > 0. The next step is to look at the D term as for 
D y^ and E = the sign of D determines whether one is dealing with fermions 
{D < 0) or bosons {D > 0) where the latter case includes (5, S), {S, V) and (V, S). 
In case oi D = only the case {S, S) remains. To further distinguish the cases {S, V) 
and (V, S) from each other, one has to consider the ratios: 



(S, V) : D/C = ^:Z., e 

(y, S) : DjB = JZlZ.. e 



-oo,-i] U [0,oo] 



0,j^] 



-2,1] 

1 22] 
2' 25 J 



0,1] 



(3.4) 



ty{ty+6)+25 

since here the dependence on the unknown couplings cancel as can be seen from 
eqs. (|A.6| ) and (|A.7| ). For these decays we have three possible ratios shown in 
Eq. (|3.4|). We see that the ratio C/B has no overlap and hence best suited ratio 
to distinguish between {S, V) and (V, S). The D/B and D/C ratios seem to be less 
useful since the intervals overlap, but except for ry = 1, e.g. my = nix, but they 
are never equal for ry 7^ 1. In the range where Ty is close to one, the SM- fermions 
become very soft and this part will be excluded because a lower cut on their energies 
is put in practice. Last but not least we note that the {S, S) can be further checked 
by the requirement that B/C = 1 as can be seen from Eq. ( [A.3|) . Therefore, indepen- 
dent of the mass ratios, one can state that all those five cases can be discriminated 
from each other. 

3.3 Impact of the third topology 

We have seen in section |2.2| that the dominant contributions stem from topologies 
1 and 2 of Table |l| in case of decays of bosons because the third topology generally 
contributes at O(e^). One might asked if one of the dimensionful couplings in the 
diagrams of the third topology can become so large to disturb the above strategy. 
Note, that in case of a new fermions higher orders in e have no impact and, thus, we 
restrict the discussion here to decays of bosons. 
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E<0 



X 



V 



>0 



k 



X 



S^ffS 4 



^S^ffV, V^ffS, S^ffS] 



X 



V^ffS 4 



■* 



v^//^^ 






■* 



//^ 






1 



Ratios: 
D/C 
C/B 
D/B 



1 



//^^ 



Figure 2: Flowchart for the strategy to discriminate different spin assignments for mf = 



using the signs of the coefficients given in section 3.1. Green (solid frame) boxes are for the 
case of taking only the leading order into account. Impact or dominance of higher order 
terms is given by the yellow (dashed frame) boxes, see text for details. 



Let us first consider the case that X is in a non-trivial SU{3) representation. Here 
only in the {S, S) there is a potentially dangerous contribution because the trilinear 
scalar coupling can in principle be of order mj. However, from Eq. ( |A.3|) we see that 
only the equality B = C gets broken, the D-term will only get a tiny contribution, 
and thus the general strategy should still work. In case of the scalar exchange in 
the third topology of the {S, V) case the coupling is momentum dependent and is of 
order mx and one is safe again. In the case X is a colour singlet all diagrams for 
the third topology in Table |I] contribute in principle. Here one has to distinguish 
two cases: (i) there are no new vector bosons or the new vector bosons do not 
belong to a new gauge group. In this case a detailed inspection of the diagrams 
shows that one arrives at the same conclusions as above, because all dimensionful 
couplings, which had not been considered before, have to be of the order mx due 
to SUii^) gauge invariance. (ii) There is a new gauge group at higher energies to 
which the intermediate vector bosons belong. In this case the SVV coupling as well 
as the masses of the vector boson will be of same order of magnitude and, thus, our 
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Table 3: 

nif = 0. 





iS,S) 


iS,V) 


iV,S) 


iV,V) 


e' 


s c 1+2 


s c/ 1+2 


s c/1+2 


s c/1+2 


B, 


+ + 


+ + 


+ + 


+ + 


c. 


± + 


± ± 


± 


± ± 


D, 


+ 00 


± 


± ± 


— — 


E, 





+ + 


+ + 


+ + 


F4 











+ + 



e^ coefficients for the case that e^, and thus also e^, are fine-tuned to vanish; 



assumption that the intermediate particle is much heavier than the decaying one 
does not hold in this case. 

3.4 Special combinations of couplings 

Up to now we have considered only the leading terms in case of decays of new bosons. 
However, it can happen that for special helicity assignments, the leading order e^ 
becomes zero. An example is the {S, S) case as can be seen in eqs. (|A.3| ) and (|A.4|) 
where the leading order is proportional to {g{r, s)n{l, s) + g{l, s)n{r, s))^ which is 
in the case of e.g. g{l, s) = n{r, s) = 0. In such cases the ]/ /m'j part of the fermion 
propagator becomes important which contributes only at 0(e'^). The question now 
is to which extent we would arrive at wrong conclusions using the strategy discussed 
so far. In table ^ we give the resulting signs of the coefficients which have to be 
compared with the O(e^) coefficients of Table 0. 

We start with the case where X is charged under SU{3). One immediately sees that 
the cases {S, S) and (V, V) are not affected. The problematic ones are {S, V) and 
(V, S) which now get a positive E as is the case of (V, V) in the leading order. Unfor- 
tunately, ratios of the other coefficients do not help if one has no further information 
on couplings and/or masses of the intermediate particles. We have marked this pos- 
sibility in Fig. ^ with the yellow boxes surrounded by red dashed lines. However, 
we want to stress that this requires fine-tuning between different couplings which 
although being unlikely, cannot be excluded on logical grounds. 
In the case that X is a SU{3) singlet the situation gets even a little bit more com- 
plicated, because now also in case of {S, S) the D coefficient is non-zero. However, 
it is still positive and, thus, it can for sure not be confused with the case of a new 
fermion. In the (V, V) case on the other hand we get in principle even more informa- 
tion as now also the F is non-zero which immediately tells us that there is a special 
combination of couplings. 

3.5 Final states containing massive SM-fermions 

Here we summarise the changes for massive SM-fermions, which in practice only is 
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iS,S) 


iS,V) 


iV,S) 


iV,V) 


e^ 


s 


c/1+2 


s 


c/ 1+2 


s 


c/1+2 


s c/1+2 


A2 








+ 


+ 


+ 


+ 


+ + 


B2 


± 


± 


± 


± 


± 


± 


± ± 


C2 


± 


± 


± 


± 


+ 


+ 


± ± 


D2 








± 


± 


+ 


+ 


± ± 


E2 




















+ + 


e" 










A, 








± 


± 


± 


± 


± ± 


B, 


± 


± 


± 


± 


± 


± 


± ± 


C3 


± 


± 


± 


± 


± 


± 


± ± 


D^ 








± 


± 


± 


± 


± ± 


E3 




















± ± 


e^ 










z. 




















— — 


A, 


+ 





+ 


+ 


+ 


+ 


+ + 


B, 


± 


+ 


± 


± 


± 


± 


± ± 


C4 


± 


+ 


± 


± 


± 


± 


± ± 


D, 


± 


— 


± 


± 


± 


± 


± ± 


E, 








± 


± 


± 


± 


± ± 


F, 




















± ± 



Table 4: Same as Table but for massive SM-fermions. 





iS,S) 


{S,V) 


iV,S) 


iV,V) 


e^ 


s 


c 1+2 


s c/1+2 


s c/1+2 


s c/1+2 


Z4 














— — 


A, 


+ 





+ + 


+ + 


+ + 


B, 


± 


+ + 


± ± 


± + 


± ± 


c. 


± 


+ + 


± ± 


± 


± ± 


D, 


+ 





± ± 


± ± 


± ± 


E, 








± + 


+ + 


± + 


F, 














+ + 



Table 5: Same as Table ^ but nif 7^ 0. 

important for top-quarks. The signs of the coefficients are given in Table ^. It turns 
out, that things hardly change but for the fact that one has to fit more coefficients. 

Comparing tables ^ and |^ one sees that the same strategy can be used in principal. 
However, for distinguishing between the (5, V) and (V, S) cases the ranges for the 
ratios of the coefficients change. Moreover, only in the ratio D/C the unknown 
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couplings cancel and we find 



(S, V) : D/C = J7^^^_, [-00, -i] U [0, 00] 



(3.5) 



where one has to use 



1 > ry + 2r/ (3.6) 

due to total energy /momentum conservation. As in the case of massless SM-fermions, 
the overlap region of the two intervals is for the case Ty — t- 1— 2rj, e.g. the kinematical 
limit, where all particles are practically at rest in the centre of mass system of X. In 
general this ratio will be either negative or much larger than 1/2 in the (S", V^ case. 
Note that we have A = for the (S", 5) case which, thus, serves as a confirmation of 
this case. Also in case of a new fermion we arrive at the same conclusions because 

sign(A4) = + , sign(54) = ± , sign(C4) = ± , sign(D4) = - , sign(E4) = 

(3.7) 

The only exception is where this decay is mediated solely by scalars in the third 
topology as in this case ^4 = 0. 

For completeness, we also give the results in the case that the leading orders vanish 
in the case of decays of bosons in Table ^ It turns out that this case is the same 
as for the case of massless SM-fermions discussed above except that now in general 
also A will be non-zero. 

3.6 Dependence on the mass of the intermediate particle(s) 

We now address the question how small e has to be so that our strategy works. 
For this we consider two examples: (I) gir) = g{l) = n{r) = n{l) = 1 where 
the leading order dominates and (II) g{r) = n{r) = 1 but g{l) = n{l) = so 
that the leading order vanishes in case of the bosonic decays and the subleading 
orders become dominant. Note, that in case of new fermions we did not manage 
to find a combination where the leading order vanishes. In all cases we have taken 
mx = 1 TeV, my = 100 GeV and nif = 0. We have checked that our results do 
not depend crucially on these values except for the cases where my gets close to 
mx which would imply soft SM fermions and experimental difficulties to observe the 
decay. 
In figures ^ and ^ we show the relative deviation 

R = -^ ^ with dFi = -— ^ , 3.8 

ai H lids 

and H denotes the limit m/ — )■ 00 and F^ the differential width for a given e. We find 
that for a decaying scalar and for a decaying fermion the deviation is always below 
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Figure 3: Ratio R for the processes X — )• ffY for scenario (i) taking mx = 1 TeV, 
nif = 0, my = 0.1 TeV. We have calculated the decay width for the following masses of 
the intermediate particle: mi = 2, 3, 5, 7, 10, 15, 30 TeV. 



20%. In case of a decaying vector particle the situation is more difficult and only 
for e < 1/5 we get i? < 0.2 for all values of s. The reason for these large deviations 
for |s| close to one is, that here the differential widths becomes zero and the rise/fall 
at the ends of the interval gets steeper the smaller e is. This also implies, that in 
the corresponding intervals for s one will observe only a few events. The situation 
improves for a decay of a vector boson if the subleading terms become dominant as 
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4: Relative deviation for differential decay width the processes X — )• ffY with 
;s (II). 



in Fig. 



In summary: our strategy should work well in all cases if e is below 1/5, in cases of 
scalar particles or spin 1/2 fermion e = 1/2 is already a reasonable value. This is for 
example a natural value for gluino decays in supersymmetric models. 



18 



4. Testing our Strategy with Monte Carlo Simulations 



We now test our strategy in Monte Carlo simulations at the parton level. For this 
one of us created data sets and the second one tried to find out the spin-assignments 
without any prior knowledge but my and {mx — rny) and their uncertainties. We 
presume that the first one is given by an independent source with a precision of 
10% and the second one with a precision of 3%. For definiteness we have taken 
mx = 1 TeV, my = 100 GeV and rrif = 0. In the following we will denote by {B, C), 
{B, C, D), {B, C, D, E) and (E, C, D, E, F) the differential width in Eq. Q where 
all but the given coefficients are zero. 



4.1 Fitting procedure 

In practice one will not have ^(s) but one will have the number of events for a given 
interval [sj, Sj + As]. For this reason we actually fit 'distributions' of the form 






-1+iAs 



l+(i-l)As 



dT ,. 
-—as 
as 



(4.1) 



where As = 2/n, n is the number of bins considered and -^ is given by Eq. (|3.1| ). 
For the creation of the 'data' we have used our model file for generic particles and 
couplings pg for the O 'Mega/ WHIZ ARD Monte Carlo generator |0|, ^ which 
contains generic particles and couplings. 

For fitting we use a linear least squares approach as described e.g. in [^ and the 
references given therein. We will exemplify this for the case of {B, C, D): 



( X\ x\ x\ \ 
X2 x\ x\ 



Vnn nn^ nn'-^ I 

•^n -^n -^n / 




X ■ 




( datai \ 
data2 

\ data„ / 



(4.2) 



where n is the number of bins. This equation can be solved by rewriting it as: 



{XfX-c= (Xfd 



(4.3) 



where c = [B, C, D)^ and d contains the 'measured' data of the differential width 
integrated over intervals of length As. There are various methods to solve this 
equation, e.g. QR-decomposition. The fit is reliable if the matrix (X)^X is well- 
conditioned, e.g. if its eigenvalues are of similar order of magnitude. After solving 
Eq. ([4.3|) we calculate the corresponding x^ 



X 



n 



E 



(Expectedj — Observedj 
Expectedj 



(4.4) 
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where n is the number of bins (= number of data points) and j the degrees of 
freedom (in this case the number of coefficients) of the fit function. We estimate the 
error for the coefficients by adding a Poisson noise to the data and determine the 
corresponding confidence interval after fitting. 

4.2 A super symmetric example 



As a first test we study a focus point scenario which is inspired by SPS2 [^: ra^ = 
3 ■ 10^ GeV, mi/2 = 3 ■ 10^ Ao = 0, tan/3 = 10 and sign(/i) > 0. We simulate for 
the gluino decay 2k and 10k events and study also the effect of different binning 
sizes, namely 10 and 50 bins. Afterwards we fit the resulting distribution to all 
possible cases. Here we have assumed that the following information on the masses 
is given: {nix ~ f^v) = 688 ± 23 GeV, my = 121 ± 12 GeV with a Gaussian 
distribution. The results are summarised in Table |^. As one does not know neither 
the absolute values of the couplings nor the intermediate masses, one has the freedom 
to normalise i? to 1 and, thus, only the other coefficients and their uncertainties are 
given. The x^ favours for both 10 and 50 bins slightly the {B, C, D) polynomial. 
However, if higher powers in s are included, one still obtains a good fit and the x^ 
will not be sufficient to discriminate between the different possibilities. This is a quite 
generic feature because usually there is a hierarchy between the non-zero coefficients: 
\B\, \C\ ^ \E\, \F\. The {B,C) case can be ruled out since B ^ C within the range 
of the error bars. The smallness of the parameters and the large errors on the E 
and F coefficients of the {B, G, D, E) and (B, C, D, E, F) models suggest that those 
values are zero. The remaining model is {B, G, D) with negative D which suggests 
(F, F). We find it encouraging that one gets already for 2000 events first information 
including that the (S", S) case can be excluded. 

4.3 Large sample tests 

In the second step we have tested our strategy for a large set of random couplings 
fixing however the kinematics to nix = 1 TeV, my = 100 GeV and all nij = 15 TeV. 
The latter number is not crucial as long as it is above 5 TeV (2 TeV) in case of 
decaying vector bosons (decaying scalars and fermions). We have generated 100 
different sets with random couplings for each of the decays {S,S), (5", V"), iV,S), 
(V, V) and (F, F) and we have generated for each set of couplings 10^, 10^ and 10^ 
events. In an ideal world one could use the strategy depicted in Fig. || without any 
problems. In reality there will be some smearing of the data from the measurement 
itself as well as from the background subtraction. To be sure that we do not miss 
anything we have slightly advanced our strategy and apply it to the following model 
hst 

{(5, S), (S, V), {V, S), (F, F), {V, V), {V, V),} (4.5) 



20 



Coef. 


x' 


(F,F)io 


x' 


(F,F)5o 


x' 


(F,F)io,2fe 


C 


5.9 


0.210 ±0.061 


1.50 


0.208 ±0.071 


1.58 


0.208 ±0.126 


C 


0.55 


0.089 ±0.077 


0.62 


0.089 ±0.076 


1.55 


0.089 ±0.155 


D 




-0.232 ±0.106 




-0.227 ±0.110 




-0.227 ±0.252 


C 


0.63 


0.079 ±0.141 


0.63 


0.067 ±0.128 


1.75 


0.067 ±0.300 


D 




-0.218 ±0.150 




-0.197 ±0.139 




-0.197 ±0.322 


E 




0.028 ±0.227 




0.060 ±0.228 




0.060 ±0.536 


C 


0.75 


0.080 ±0.149 


0.64 


0.077 ±0.160 


1.14 


0.077 ±0.401 


D 




-0.215 ±0.456 




-0.169 ±0.340 




-0.169 ±1.108 


E 




0.025 ±0.260 




0.035 ±0.311 




0.035 ± 0.738 


F 




-0.006 ±0.516 




-0.046 ± 0.443 




-0.046 ±1.315 



Table 6: Testing of the SUSY focus point taking 10 bins in case of 2000 events and 
10 and 50 bins for 10000 events. Input for the fit is: {mx — fny) = 688 ± 23 GeV, 
my = 121 ± 12 GeV with gaussian distribution, mi = 3026 GeV. The coefficients are 
normalised such that B = 1 and the uncertainties are at 3a. The analytic values for the 
coefficients aie B = I, G = 0.123, D = -0.188. 

and the corresponding differential widths. {V, V)4 denotes the case, where only the 
4th order in e of (V,V) remains. We have tested for the following criterions: 

1. B^C ^ remove {S, S) 

2. D>0 ^ remove {F,F); D <0^ remove (S, V), {V, S) 

3. C/B in (S,V) interval — > if not remove {S, V) 

4. C/B in (V,S) interval — )■ if not remove (V", S) 

5. D/B in (S,V) interval — )■ if not remove {S, V) 

6. D/B in (V,S) interval — ?■ if not remove (V, S) 

7. D/C in (S,V) interval — > if not remove {S, V) 

8. D/C in (V,S) interval — ?■ if not remove (V", S) 

9. if E < in {V, V)/{V, V)4, ^ remove {V, V)/{V, V)^ 

10. if F < in (V, V)^, -> remove (V, V)^ 

11. Optional: Remove all models with x^ > 3 

12. Optional: Remove {V, V)/{V, V)^ ii E,F < 0.001 respectively 
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Every time one criterion could be applied or not fulfilled within the range of the error 
bars, the corresponding model was cancelled from the list in Eq. ( [4.5|) . 
In table |^ we have summarised our results for these different Monte Carlo data sets 
where we give the number of the remaining models after going through the different 
criteria. The first row e.g. means, that we started with 100 different data sets for 
the {S, S) decay, applied our tests and after that, 100 (5, 5*) models remained, of 
{S, V), 20 of (V, S), 100 of {F, F) and so on. There were no data sets where only a 
wrong model remained. 





(S,S) 


(S,V) 


(V,S) 


(F,F) 


(V,V) 


(V,V4) 


10"^ events: 


(S,S): 


100 





20 


100 


100 


100 


(S,V): 





99 





1 


100 


99 


(V,S): 








99 


99 


100 


98 


(F,F): 





5 





100 


99 


99 


(V,V): 





66 





78 


100 


100 


10^ events: 


(S,S): 


100 








100 


100 


100 


(S,V): 





99 








100 


100 


(V,S): 








100 


100 


100 


100 


(F,F): 











100 


100 


100 


(V,V): 





66 





78 


100 


100 


10^ events: 


(S,S): 


100 








100 


100 


100 


(S,V): 





98 








99 


100 


(V,S): 








100 


100 


100 


100 


(F,F): 











100 


100 


100 


(V,V): 





10 





61 


100 


100 



Table 7: Number for the remaining models for 100 input models each where the various 
criterions are applied using 3a uncertainties on the coefficients. The masses are chosen 
as my = 100 lb 10 GeV and the mass difference {mx — rny) = 900 it 30 GeV . The bold 
numbers are the correct model fits. 

The obtained results can be understood as follows: (i) It is easier to fit a polynomial 
which has low powers of s by a higher order polynomial if there is smearing than vice 
versa, (ii) The number of criterions depend on the decay mode, e.g. it is easier to 
exclude (5*, V) where 3 criterions are at hand than {V, V) where only one exists, (iii) 
The modulus of the coefficients E and F is usually up to two orders of magnitude 
smaller than the modulus of the other coefficients but the absolute uncertainty is 
roughly the same for all coefficients. In particular the third item implies that it will 
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(S,S) 


(S,V) 


(V,S) 


(F,F) 


(V,V) 


(V,V4) 


10^ events: 


(S,S): 


99/100 


0/0 


19/20 


99/100 


99/85 


97/90 


(S,V): 


0/0 


96/96 


0/0 


0/1 


93/95 


90/89 


(V,S): 


0/0 


0/0 


97/97 


97/99 


99/88 


96/81 


(F,F): 


0/0 


5/5 


/OO 


98/100 


98/89 


97/87 


(V,V): 


0/0 


46/66 


0/0 


31/78 


95/95 


92/90 


10^ events: 


(S,S): 


100/99 


0/0 


0/0 


100/100 


97/63 


98/63 


(S,V): 


0/0 


99/95 


0/0 


0/0 


97/81 


94/72 


(V,S): 


0/0 


0/0 


100/100 


100/100 


100/66 


99/64 


(F,F): 


0/0 


0/0 


0/0 


99/100 


98/66 


98/61 


(V,V): 


0/0 


46/66 


0/0 


31/78 


95/95 


92/90 


10^ events: 


(S,S): 


100/100 


0/0 


0/0 


100/99 


100/14 


96/23 


(S,V): 


0/0 


95/98 


0/0 


0/0 


93/54 


90/33 


(V,S): 


0/0 


0/0 


99/100 


99/100 


98/17 


98/28 


(F,F): 


0/0 


0/0 


0/0 


100/100 


98/22 


97/25 


(V,V): 


0/0 


0/10 


0/0 


0/61 


93/100 


97/64 



Table 8: Same as Table but taking into account either the optional criterion 11 ( exclude 
X^ > 3) or criterion 12 (small E,F) corresponding to the first and second number given at 
the various entries. 



be rather difficult to exclude a positive E and F in practice if only one decay channel 
is considered. 

Additionally we have also looked at the optional criteria item |Tl| and |12| separately 
and give the resulting numbers in Table |^. The ffist number of 99/100 (as e.g. in the 
{S, S) case) denotes the number remaining after applying the x^ criterion, the second 
number the same with the small E,F criterion. The results are (i) the x^ > 3 criterion 
is most useful if the underlying decay is (V, V)/(y, V")4, since the polynomials with 
a lower order have a large x^- This is e.g. reflected in the {V,V) decay with 10k 
events, and the fitted {S, V) and {F, F) polynomials. Here the number of remaining 
processes are reduced by applying the x^ test from 66 {{S, V), Table |^ to 46 (Table |^) 
and for {F,F) from 78 to 31. (ii) The criterion for small E,F works very good for 
high statistics since the ffi gives values close to 0. This is reflected in the last row 
of Table P, where the underlying process is {F,F). After applying of the common 
criterions, 100 of (V, V) and 100 of {V, V)^ models remain (see Table |^). But after 
applying additionally the criterion item |T2| this number is reduced to 22 (V, V) and 
25 (V, V)^. However, the same argument does not remove any of the (V, V) models, 
if (V", V) is the underlying process. For a smaller number of events, there are only 
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a view coefficients E,F smaller then 0.001 so only a smaller number of models are 

removed as e.g. in Table ^ for 10k events, 88/81 remaining (V, V")/(V, 1^)4 models 

for the underlying (V, S) process. 

Since our decision making strategy depends upon the correct error estimation of 

the coefficients, we have independently checked that the Markov-Chain-Monte-Carlo 

(MCMC) method for the data fitting which yields roughly the same errors on the 

coefficients. 

5. Conclusions 

In this paper we have investigated three-body decays of the form X — )■ f fY where 
X and Y are new particles and / are SM-fermions with the aim to determine the 
spins of the new particles. Here we assumed that F is a DM candidate and escapes 
detection in a typical LHC detector. 

We have studied in detail the differential width as a function of the invariant mass 
of the SM-fermions for the case of heavy intermediate particles with mass mi and 
expanded the width in powers of the ratio e = mx/rni. It turns out that general 
properties such as signs or various ratios of the resulting coefficients depend on the 
spin assignments of X and Y. From this we have developed a strategy for the spin 
identification discussing various cases and testing it on large samples of arbitrary 
coupling assignments. Here it turns out that one is able to exclude several spin 
assignments but one does not get necessarily find a unique solution once one has to 
deal with noisy data. 

Although we did not find a unique solution we are convinced that the proposed 
method will be useful in practice for the following reasons: (i) We have only in- 
vestigated one particular decay channel. However, in general several channels will 
be open which can be combined, (ii) In the same spirit: we have only investigated 
one decaying particle. In practice, e.g. if supersymmetry or extra dimensions are 
realized in nature, several distinct new particles will be produced which eventually 
have to decay into the lightest of the new ones if a parity like i?-parity or KK-penity 
is realized. Therefore, one will have several different possibilities to determine the 
spin of the invisible particle Y. (iii) Our information can be combined with other 
observables, e.g. with cross section information. However, here one most likely will 
have to assume a certain representation to which a particular new particle belongs, 
e.g. if it is an SU{3) triplet or octet. 
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A. Analytic Results for the Coefficients 

The coefficients of the differential decay widths are given below. We restrict ourselves 
to the case of massless SM-fermions, implying that Z = A = and, thus, the 
differential widths read as 

^ = T^^F^I {B + CS + Df + Es' + Fs') (A.l) 

as [Znj'^ 256 nix 

A.l Decays of new bosons 

The coefficients are shown with all possible diagrams and vertices in Table |l|. For the 
definition of e and the various Tj see Eq. (p.2|). We give the various orders separately, 
e.g. 

4 

B = Y,Bje^ (A.2) 

i=2 
For brevity, we only explicitely write out the higher orders for [S, S). 

S ^ f7S: 

B2 = 128e\TY - Ifigir, s)n{l, s) + g{l, s)n{r, s)f 
C2 = l28e'{TY - Ifigir, s)n{l, s) + (?(/, s)nir, s)f 
D2 = (A.3) 

fig = 64eVc(ry - iyc{.s){g{r, s)n{l, s) + g{l, s)n{r, s)){s{l) + s{r)) 
Cs = 64eVc(ry - l)'cis){g{r, s)n(/, s) + g{l, s)n{r, s)){s{l) + s{r)) 
Ds = (A.4) 

B4 = y e'(T-y - If {Ugir, sfn{l, sf{TY + 1)' + I2g{l, sfn{r, sf{TY + 1)' 
+24^(/, s)g{r, s)n{l, s)n{r, s)(ry + 1)^ + 3T^c{sfs{lf + 3r^c{sfs{rf 
+4c{v)M^y + QTYc{v)h{iy + c{vfv{lf + 4c{vfv{ry 
+6TYc{vyv{rf + c{vfv{ry) 
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s 



C4 = -y e^(^K - 1)' (-48ry(7(r, s^nil, s)^ - 9QTYgil, s)gir, s)n{r, s)n(/ 
-48ry(7(/, s)2n(r, s)^ - 3T^c{sfs{lf - 3T^c{sys{rf + 24c{v)Miy 
+4TYc{vyv{l? + 2c{v)Ml? + 24c{vyv{rf 
+4ryc(t;)\(r)2 + 2c{vyv{r)'^) 

16 

D, = -y6'(ry - 1)' (- {vilf + i;(r)^) civ)' + I2g{r, sfn{l. sf 

+12g{l, sfn{r, sf + 245f(/, s)g{r, s)n{l, s)n{r, s)) (A.5) 

Moreover, we get Ej = in all orders considered. 
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S ^ ffV: 

64 
^2 = ^ {gir, sfn{l, vf + g(l, sfn{r, vf) e^ry - if {2hT^ + 6ry + l) 

1 28 
C2 = ^ (^(r, s)2n(/, vf + (?(/, s)2n(r, t;)^) t^ry - if (llr^ - 2ry - I) 

fi4 
^2 = ^ {9{r. sfnil, vf + g{l, sfn{r, vf) e\ry - if (A.6) 

V ^ flS: 

64 
^2 = y (^(r, t;)=^n(/, s)^ + (?(/, t;)2n(r, s)^) e2(^^ _ if (^2 ^ g^^ ^ 25) 

C, = -^ (^(r, vfn{l, sf + (?(/, t;)2n(r, sf) e\ry - if {4 + 2ry - ll) 

64 
D2 = \ (^(r, vfn{l, sf + (?(/, t;)2n(r, s)^) ^^(ry - 1)^ (A.7) 

32 

^2 = ^-^e2(ry - 1)2 [g{r, vf (3r^ + IQt^ + 54r^ + 16ry + 3) n(/, t;)^ 

-g{l, v)g{r, v)n{r, v) (3r^ + 20^ - Sr^ + 20ry + 3) n(/, t;) 

+g{l, vfn{r, vf (3r^ + 16t^ + 54r^ + 16ry + 3)) 

32 

C2 = -^^\ry - If {g{T, vf (5r^ + 4r^ - 46^ + 4ry + 5) n(/, vf 

-g{l, v)g{r, v)n{r, v) (7r^ + 20^ + 34^ + 20ry + 7) n{l, v) 
+g{l, vfn{r, vf (5r^ + At^ - 46r2 + 4ry + 5)) 

32 

D2 = ^t\ry - 1)^ [g{r, vf (r^ - 6ry + l) n{l, vf 

—g{l, v)g{r, v)n{r, v) (pTy + 6ry + 5) n(/, v) 
+g{l,vfn{r,vf {4 - 6ty + I)) 

32 

^2 = ^ (fl'(^, vfn{l, vf + g{l, v)g{r, v)n{r, v)n{l, v) + g{l, vfn{r, vf) e'^{Ty - if 

(A.8) 

A. 2 Decays of new fermions 

As noted before, in this case only the 4th order in e contributes. We spht the 
various coefficients according to the different topologies considered, e.g. the scalar 
contributions to B4 are 

54 = 5, + 5; + B.,,m (A.9) 

where 

Bg : top. 1+2 with intermediate scalars 



26 



B'g : top. 3 with intermediate scalars 

Bgfn '■ interference term of top. (l+2)+3 with intermediate scalars 

For intermediate vector bosons the index v is used. For the interference terms be- 
tween scalars and vector bosons the index (sv) is used in case of topologies 1+2, 
(svl) for the scalars of topology 1+2 and vector bosons of topology 3, (slv) for the 
scalars of topology 3 and vector bosons of topology 1+2. Moreover we find that the 
interference vanishes if both, scalars and vector bosons stem from the third topology 
because we have ruf = for the SM-fermions. 

Intermediate scalars: 



64 
Bs = y (^^^ - 1)' (2 {(Mry + 6) + l)n(/, sf + {ty{ty + 3) + l)n(r, sf) g{l, sf 

-{TYiTY + 6) + l)g{r, s)n{l, s)n{r, .s)g{l, s) + 2g{r, s)^ {{TYirY + 3) + l)n(/, sf 
+ (ry(ry + 6) + l)n(r,s)2)) 

64 
Cs = -y (ry - 1)' (-2(?(/, s)(?(r, s)n(/, s)nir, s)irY + 1)^ 

+g{r, sf (((ry - 4)ry + l)n{l, sf + {{ty - 10)ry + l)n{r, s^) 

+gil, sf {{{ty - 10)ry + l)n(/, sf + {{ty - A)ty + l)n(r, sf)) 

64 
Ds = --{ty - If {{nil, sf + n(r, sf) g{l, sf 

+g{r-, s)n{l, s)n{r, s)g{l, s) + g{r, sf {n{l, sf + n{r, sf)) (A. 10) 



B'^ = 32{TY-lf{{TY + lfd{l,sf 

+8ryrf(r, s)d{l, s) + {ty + l)'rf(r, sf) {s{lf + s{rf) 
Ci = 128(ry - 1) Vy (rf(/, s) + rf(r, s)f {silf + sirf) 
D'^ = -32(ry - 1)^ (rf(/, s)2 + rf(r, sf) {silf + s(r)2) (A.ll) 

Bs,„, = -32(ry - l)2((7(r, s)n{l, s) + (?(/, s)n(r, s)) {d{r, s) {s{1){ty + 1)' + 4TYs{r)) 

+rf(/,s)(s(r)(ry + l)2 + 4rys(/))) 
C,,„ = -128(ry - lfTY{g{r, s)n{l, s) + g{l, s)n{r, s)){d{l, s) + d{r, s)){s{l) + s(r)) 
D,,„ = 32(ry - lf{g{r, s)n(/, s) + <?(/, s)n{r, s))(rf(r, s)s(/) + rf(/, s)s{r)) (A.12) 

Intermediate vector bosons: 

256 

S„ = — (T-y - 1)^ {Qg{l,v)g{r,v)n{l,v)n{r,v){TY + lf 

+g{l, vf {2iTYiTY + 6) + l)n(/, t;)^ + 3(ry + lfn{r, vf) 
+gir, vf (3(ry + if nil, vf + 2{TYiTY + 6) + l)n(r, t;)^)) 
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a = —{ry - If {{l2Tyn{r,vf - {{ry - 10)ry + l)n{lvf) g{lvf 

+2ATYgir, v)n{l, v)n{r, v)g{l, v) + g{r,vf (l2rYn(/,w)^ - {{ry - 10)ry + l)n{r,vf)) 

D, = -^{ry - 1)^ ( {nil, vf + 3n(r, vf) g{l, vf 

+6(7(r, v)n{l, v)n{r, v)g{l, v) + g{r, vf (3n(/, vf + n{r, vf)) (A. 13) 

2S6 

K = —i^y- 1)' HMry + 3) + l)d{l, vY - 6ryd{r, v)d{l, v) 

HMry + 3) + l)d{r, vf) {v{lf + v{rf) 
C. = -^(ry - 1)^ {{{ry - 4)ry + 1^/, vf 

+ l2Tyd{r, v)d{l, v) + {{Ty - A)Ty + l)d{r, vf) {v{lf + v{rf) 
D'. = -^{ry - If {d{l.vf + d{r,vf) {v{lf + v{rf) (A.14) 



512 
B,,m = —{ry - lf{ry{Ty + 6) + l){d{l, v) - d{r, v)) 

{g{l, v)n{l, v)v{l) — g{r, v)n{r, v)v{r)) 
C,,„^ = - — {ry - lf{{ry - 10)ry + l){d{l,v) - d{r,v)) 

{g{l, v)n{l, v)v{l) — g{r, v)n{r, v)v{r)) 

256 
Dv,m = —{ry - lf{d{l, v) - d{r, v)){g{l, v)n{l, v)v{l) - g{r, v)n{r, v)v{r)) 

(A.15) 

Interference terms between scalars and vector bosons: 

512 
Bsv = :^{ry - lf{n{l,v){3Tyg{l,s)g{r,v)n{r,s) + g{r,s){3Tyg{r,v)n{l,s) 

+ {ry{ry + 6) + l)g{l, v)n{r, s))) + {{ry{Ty + 6) + l)g{l, s)g{r, v)n{l, s) 
+ 3ry5f(/, v){g{r, s)n{l, s) + g{l, s)n{r, s)))n{r, v)) 

256 
Csv = —{ry - lf{n{l,v){g{r,s){{{Ty - 10)ry + \)g{l,v)n{r,s) 

- Qryg{r, v)n{l, s)) - QTyg{l, s)g{r, v)n{r, s)) + {{{ry - 10)rY + \)g{l, s)g{r, v)n{l, s) 

- QTyg{l, v){g{r, s)n{l, s) + g{l, s)n{r, s)))n{r, v)) 

256 

Dsv = -ir{'^y ^ l)''(fl'(^> v)g{r, s)n{l, v)n{r, s) + g{l, s)g{r, v)n{l, s)n{r, v)) (A.16) 



256 
Bs.i = - — {ry - lf{ry{ry + 6) + l){d{l, v) - d{r, v)) 

{g{r, s)n{r, s)v{l) - g{l, s)n{l, s)v{r)) 
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Csvi = ^{ty - Ifiiry - 10)ry + l){d{l,v) - d{r,v)) 

{g{r, s)n{r, s)v{l) - g{l, s)n{l, s)v{r)) 
Ds.i = ^^{ry - l)\d{Uv)~d{r,v)){g{r,s)n{r,s)v{l)-g{l,s)n{l,s)v{r)) (A.17) 



5,i„ = 128(ry - lf{g{r, v)n{l, v) + g{l, v)n{r, v)) {d{l, s) (s(/)(ry + 1)^ + 4rys(r)) 

+rf(r,s)(s(r)(ry + l)2 + 4rys(/))) 
Csiv = 512(ry - 1) Vy (5((r, v)n{l, v) + g{l, v)n{r, v)) 

(rf(/,s) + rf(r,s))(s(/) + s(r)) 

D,i„ = -128(ry - l)\g{r, v)n{l, v) + g{l, v)n{r, v)){d{l, s)s{l) + d{r, s)s(r)) 

(A.18) 
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